Low-Redshift Cosmic Baryon Fluid on Large Scales and She-Leveque Universal Scaling 
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We investigate the statistical properties of cosmic baryon fluid in the nonlinear regime, which is 
crucial for understanding the large-scale structure formation of the universe. With the hydrodynamic 
simulation sample of the Universe in the cold dark matter model with a cosmological constant, we 
show that the intermittency of the velocity field of cosmic baryon fluid at redshift z — in the scale 
range from the Jeans length to about 16 h _1 Mpc can be extremely well described by She-Leveque's 
universal scaling formula. The baryon fluid also possesses the following features: (1) for volume 
weight statistics, the dissipative structures are dominated by sheets, and (2) the relation between 
the intensities of fluctuations is hierarchical. These results imply that the evolution of highly evolved 
cosmic baryon fluid is similar to a fully developed turbulence. 

PACS numbers: 98.65.Dx, 47.27.Gs 
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Introduction. — The formation and evolution of large- 
scale structure of the universe is governed by the gravi- 
tational clustering of cosmic matter, in which about 72% 
is dark energy, 24% cold dark matter, and 4% baryon. 
Most baryonic matter is in the form of gas. Therefore, 
the evolution of the density and velocity fluctuations of 
cosmic baryon fluid is dominated by the underlying grav- 
itational potential of dark matter. In the linear regime, 
the baryon fluid follows the mass and velocity fields of 
collisionless dark matter point by point. In the nonlinear 
regime, however, as first pointed out by Shandarin and 
Zeldovich in their early study of structure formation, the 
dynamical behavior of cosmic matter clustering on large 
scales is similar to turbulence 1] . They emphasized that 
the motion of self-gravitation matter in the expanding 
universe is like that of noninteracting matter moving by 
inertia. In other words, cosmic matter underwent a scale- 
free evolution somewhat like fully developed turbulence 
in inertial range. 

Nevertheless, the dynamical difference between the 
structure formation of cosmic matter and turbulent flow 
in incompressible fluid is obvious. The latter is rotational 
in general 0, while the former is irrotationaL because 
vorticities do not grow during the clustering In tur- 
bulence, energy passes from large to the smallest eddies, 
and finally dissipates into thermal motion, while cosmic 
baryonic gas falls into massive halos to form structures, 
including light-emitting objects. 

Yet, the turbulencelike behavior of cosmic baryon fluid 
has been gradually noticed in the last decade. Although 
the evolution of cosmic baryon fluid is governed by the 
Naiver-Stokes equation, the dynamics of growth modes 
of the fluid can be sketched by a stochastic force driven 
Burgers' turbulence 3 
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scale function of the cosmological expansion, and v is 
from the Jeans diffusion. The stochastic term <f> is actu- 
ally the gravitational potential of the dark matter mass 
density perturbation, i.e., V 2 ^ = 47rGa 2 pdm<5dm, and 
5dm = (Pdm - Pdm)/Pdm- The stochastic property of <5 dm 
is given by the correlator (#dm(k)<5dm(k')) = -P(k)<5t,k' 
and the power spectrum P(k) oc k~ a . For the flat 
cold dark matter model with a cosmological constant 
(ACDM), the index is initially a — 2.3 in the scale range 
considered. A basic feature of the fluid described by 
Eq.(l) is that the Burgers' turbulence is developed when 
the Reynolds number is large For instance, the prob- 
ability distribution functions of the velocity differences 
across a distance r, Sv r = {[v(x + r) — v(x)] • r/r}, are 
scaling, when r is larger than the Jeans length These 
results indicate that some dynamical features of the cos- 
mic baryon fluid at low redshifts in the nonlinear regime 
are turbulencelike. 

Therefore, an important question is how well and why 
the dynamical state of cosmic baryon fluid at low red- 
shifts can be described as turbulence. Although the cos- 
mic baryon fluid has been a central topic in cosmology 
for a long time Q, no turbulencelike behavior has been 
studied. In this Letter, we investigate this problem in the 
context of the universal scaling law of fully developed tur- 
bulence. Using numerical samples, we first show that the 
universal scaling does exist in the cosmic baryon fluid, 
and then analyze the underlying physics of this scaling. 

Intermittent exponent. — The universal property of 
fully developed turbulence is measured by the structure 
functions S p (r) and intermittent exponent Cp, defined as 



S p (r) = (6vP) ~ r 



(2) 



where tp is velocity potential, v 



(1) 

(l/a)V<p, a is the 



Based on dimensional argument of hierarchical evolution, 
Kolmogorov in 1941 predicted that for fully developed 
turbulence on scales of inertial range, the intermittent 
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exponent is £ p = p/3 Experimental and numeri- 
cal results do not, however, support the p/3 law. It 
should be attributed to intermittency; i.e., turbulence 
field in the inertial range is characterized by stronger 
non-Gaussianity on smaller scales js). A remarkable de- 
velopment was made by She and Leveque 0] (SL here- 
after). They proposed that the non-Gaussian behavior of 
fully developed turbulence is determined by the hierarchi- 
cal structure originated from the Navier-Stokes equation, 
and the p/3 law should be replaced by 

C P /Ca = [1 - C(l - (3 3 )} P /3 + C(l - (3) 

where C is the Hausdorff codimension of the most dissi- 
pative structures, and parameter j3 is given by the hierar- 
chical evolution (see below). The SL formula Eq.(3) is in 
excellent agreement with various experiments of turbu- 
lence 1 1 1| , including also turbulence in compressible fluid 
[l2T ]. The SL scaling law of structure function is consid- 
ered to be universal for characterizing the fully developed 
turbulence. 

To investigate the universal scaling of cosmic baryon 
fluid, we use the cosmological hydrodynamic simulation 
samples produced by the code WIGEON (Weno for Inter- 
galactic medium and Galaxy Evolution and formatiON) 
131. This is a hybrid hydrodynamic/ TV-body simula- 
tion, consisting of the WENO algorithm 14] for bary- 
onic fluid, and ./V-body simulation for particles of dark 
matter. The baryon fluid obeys the Navier-Stokes equa- 
tion, and is gravitationally coupled with collisionless dark 
matter. We use the ACDM cosmological model with pa- 
rameters given by the recent observations of the cosmic 
microwave background radiation |15|. The linear power 
spectrum of mass density perturbations is taken from the 
fitting formulas of Eisenstein & Hu ^(| . The atomic pro- 
cesses in the plasma of hydrogen and helium of primordial 
composition, including ionization, radiative cooling and 
heating, are modeled in the same way as in 0. 

The simulations were performed in a periodic, cubic 
box of size 64 ft, _1 Mpc with a 512 3 grid and an equal 
number of dark matter particles. The simulations start at 
a redshift z = 99. A uniform UV-background of ionizing 
photons is switched on at z = 6 to heat the gas and 
reionize the universe. The temperature of the baryonic 
gas generally lies in the range 10 4 — 10 6 K, and the speed 
of sound v s in the baryonic gas is only a few km s _1 to 
a few tens km s . The Jeans length Aj yields a term 
like the viscosity \i ~ v s Xj. On the other hand, the bulk 
velocity of the baryonic gas is of the order of hundreds 
of km s _1 0]. Therefore, the Reynolds number would 
be larger than ~100 if the scales under consideration are 
larger than the Jeans length Aj, which is in the range 
- 0.1 - 0.3 h" 1 Mpc for redshifts z < 4 [l||. 

These samples are successful to model the power spec- 
trum and non-Gaussian features of the observed trans- 
mitted flux of Lya 0]. Hence it would be suitable 
to study whether the universal scaling law Eq.(3) is 
available for cosmic baryon fluid. For this purpose, we 
randomly sampled 10,000 one-dimensional sub-samples, 
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FIG. 1: Intermittent exponent ( p of cosmic baryon fluid at 
redshift z — 0. The means and error bars are from the average 
over samples on scale range from r — 1 to 16 h -1 Mpc. The 
solid line is the She-Leveque formula Eq.(2) with C — 1 and 
P 3 = 1/3. 

with each one containing 512 data points. At each point 
of a line, the peculiar velocity and mass density of the 
baryonic gas is recorded. 

We calculated the moments of velocity difference 6v r 
for r from 1 to 16 h _1 Mpc, of which the lower limit 
is larger than the Jeans length, and the upper limit is 
due to the size of the simulation box. The intermittent 
exponent £ p for the sample at redshift z — is shown 
in Fig.l. The error bars are the variance of the samples. 
The intermittency of cosmic baryon fluid is excellently 
fitted by the scaling Eq.(3) on all orders considered, if 
the parameters of SL formula [Eq.(3)] are taken to be 
C = 1 and /3 3 = 1/3. 

Our sample is v ery different from experimental and 
numerical samples 20] used to test the SL formula. The 
success shown in Figure 1 further supports that the hy- 
potheses used to derive the SL formula Eq. (3) are univer- 
sal so that it would also be available for cosmic baryon 
fluid. In this sense, the non-Gaussian behavior of the 
baryon fluid in an expanding universe at redshift z = 
can essentially be described as a fluid satisfying the uni- 
versal scaling for fully developed turbulence. 

Singular dissipative structures. — Simply speaking, the 
underlying physical picture of the SL formula is as fol- 
lows. In the inertial range, the fluid evolution is governed 
by scale-covariant interactions from the Navier-Stokes 
equations. The kinetic energy of fluid is dissipated in 
singular structures, and the energy dissipation on differ- 
ent scales satisfies hierarchical relation. 

Although the clustering of baryon fluid is governed by 
the gravity of the background mass field, the expansion 
of the Universe eliminates the gravity of the uniformly 
distributed matter. The peculiar motion of baryon fluid 
feels only the gravity given by the random fluctuations 
of the distributed dark matter. In a nonlinear regime, 
baryon fluid decouples from dark matter. It is similar to 
the decoupling of a passive substance from the underly- 
ing field during nonlinear evolution [2l| . The decoupling 



FIG. 2: Left: V • v contour of baryonic fluid for a slice with 
area 16 x 16 h -2 Mpc 2 and thickness 2 h _1 Mpc at z = 0. 
Right: the field of Vdm — v for the same sample as the left 
one. The background is the mass density contrast of baryon 
fluid log 5. 

leads to the velocity of baryon fluid generally being less 
or not larger than the velocity of dark matter Q . There- 
fore, baryon fluid is actually driven by a random force 
of the dark matter gravity. The potential of the random 
gravity is of power law, i.e., scale free. In the nonlinear 
regime, the power of fluctuations transfers from larger 
to smaller scales. Therefore, in the scale range of 1 to 
about 10 h _1 Mpc the evolution of cosmic fluid is similar 
to turbulence in inertial range. 

Since C = 1 gives a best fitting (Fig.l), the singular 
dissipative structures should be two-dimensional sheets. 
It is well known that in structure formation sheets are 
dominant in terms of volume weight statistics. For in- 
compressible fluid, the singular structures are traced by 
the amplitude of vorticity |V x v| 22 1. Since the Burgers 
turbulence is caused by shock wave Q , the singular dissi- 
pative structures are traced by the divergence of baryon 
fluid V • v. Figure 2 presents a typical V • v contour of 
baryonic gas in a slice at z = 0. One can clearly see the 
filamentary structures in two dimensions, and therefore, 
they are most likely to be sheets in three dimensions. 
In Fig. 2, we also plot the difference of velocity fields be- 
tween dark matter and baryon matter. It shows that the 
decoupling between the velocity fields of baryon fluid and 
dark matter occurred on the entire field. 

Hierarchical relation. — The parameter (3 of Eq. (3) de- 
scribes the hierarchical evolution of the cosmic baryon 
matter. The hierarchical clustering is well known in the 
theory of large-scale structure formation of the Universe. 
It has been proposed that the hierarchical clustering is 
given by the relation (Sp^) oc (Spf) 1 * -1 , where dpi is the 
density fluctuation field of the cosmic matter on scale 
I |2^| . In hydrodynamics, this relation is the so-called 
pure hierarchical relation j24j. It has been shown, how- 
ever, that the observed intermittency is not subject to 
the pure hierarchical relation if the proportional coeffi- 
cient between (Spf) and (Spf) 1 ^ -1 is scale independent 

M 

In the universal scaling scheme, the hierarchy is not 
described by mass density fluctuation 5p, but by a hier- 
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FIG. 3: The log- log relation between F p+ i(r)/ -FbfV) and 
F p (r)/Fi(r) for the cosmic baryon fluid sample. The slope 
is /?' = 0.69. 

archical relation of structure functions as 

F p+1 (r) = A p F p (rf F^r) 1 ^' , (4) 

where F p (r) = S p +i(r)/ S p (r), and number A p are scale 
r independent. F p+ i(r) describes the intensity of fluctu- 
ations. The larger the p of F p (r), the higher the intensity 
of fluctuations. The most intermittent structures are de- 
scribed by Fqo (r) . Equation (4) describes the hierarchical 
relation between fluctuations with different intensities. 
Equation (4) is invariant with respect to a translation 
in p. The parameter /?', which should be in the range 
< (3' < 1, is to measure the degree of intermittency 
of turbulent flow. The smaller the /3', the stronger the 
intermittency. For f3' = 1, the field is not intermittent. 
From Eq.(4), we have 

F p+1 (r) _ A p fF p (r)y' 
F 2 (r) Ax UiW/ ' 

If A p is p independent, we should have 

In F p+1 (r)/F 2 (r) = j3' \n F p {r) / F^r) + const. (6) 

Substituting Eq.(3) into Eq.(6), we have (3' = (3. Equa- 
tion (6) does not contain term F^r), and therefore, it 
can be directly used to test the assumption of hierarchy. 
Figure 3 presents the relation of \og[F p+ i(r) / F 2 (r)] vs 
log[F p (r) / Fi(r)], which is a perfectly straight line. The 
slope is (3' — 0.69 ± 0.02, which is also perfectly in agree- 
ment with the parameter (3 = (1/3) 1 / 3 = 0.69. There- 
fore, the assumption of a constant A p (p independence) 
is tenable. The hierarchical relation Eq.(4) with constant 
A p and (3 = 0.69 provides a good description of the hi- 
erarchical structures of cosmic baryon fluid on the scales 
of the inertial range. 

Discussion and conclusion. — We showed that all the 
moments of the fluctuations of velocity field of cosmic 
baryon fluid at z = on scales larger than the Jeans 
length have a power-law dependence on scales, and the 
intermittent exponents obey the universal law, which de- 
pends only on (1) the dimension of the most singular 
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dissipative structures and (2) the hierarchical relation be- 
tween structures with various intensities of fluctuations. 
This result strongly indicates that the highly nonlinear 
evolution would lead to the cosmic baryon fluid reach- 
ing a statistically quasicquilibrium state satisfying the 
universal scaling as that of fully developed turbulence. 
It should be emphasized that the turbulence described 
by Eq.(l) is stirred at all scales, not only at the largest 
scale as in the conventional model of turbulence. The SL 
scaling is reasonable not only for the conventional energy 
cascade but also the case of Eq.(l). 

In view of this picture, we can say that in the highly 
nonlinear regime, the statistical properties, especially the 
intermittent behavior, of the velocity fluctuations are ac- 
tually independent of the details of the dissipative pro- 
cesses. The state depends only on the dimension of dis- 
sipative structures and the hierarchical relation index. 



We believe that some statistical features, which have al- 
ready been recognized in large-scale structures, would be 
directly in consequence of the universal properties of cos- 
mic baryon fluid. 

The velocity and mass density field on scales larger 
than the Jeans length are the basic environment of the 
formation of the luminous objects. Therefore, the univer- 
sal properties of the nonlinear regime provide common 
frame of studying the dynamical and thermodynamical 
evolution of structure formation. For instance, the red- 
shift dependence of the scaling behavior would be useful 
for understanding the evolution of cosmic clustering. The 
details will be reported in the near future. 
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